Spin polarized tunneling in ferromagnet/unconventional superconductor junctions 
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We study tunneling in ferromagnet/unconventional superconductor (F/S) junctions. We include the 
effects of spin polarization, interfacial resistance, and Fermi wavevector mismatch (FWM) between 
the F and S regions. Andreev reflection (AR) at the F/S interface, governing tunneling at low bias 
voltage, is strongly modified by these parameters. The conductance exhibits a wide variety of novel 
features as a function of applied voltage. 
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Spin polarized transport and tunneling between ferro- 
magnetic and superconducting materials has become a 
vigorously pursued area of research. The studies per- 
formed [Q-pJ have implications for the understanding of 
unconventional superconductivty and for the develop- 
ment of devices using spin polarized current (TJ|^||. A 
recent study of a ferromagnet /unconventional supercon- 
ductor (F/S) system has revealed [El a differential con- 
ductance dip at zero bias (ZBCD), attributed to suppres- 
sion of Andreev reflection (AR) as a consequence of high 
spin polarization in the ferromagnet. While there is sub- 
stantial work Jl],||jl]|| on the interplay of ferromagnetism 
and superconductivty in tunneling properties involving 
s-wave superconductors, there is still no adequate theory 
for these phenomena in the context of unconventional su- 
perconductivity. Replacing an s-wave superconductor by 
an unconventional one in a normal metal/superconductor 
(N/S) structure can drastically alter the conductance 
spectrum and important changes should also oc- 

cur when such replacement is made in a F/S structure. 
For N/S junctions, zero bias conductance peaks (ZBCP) 
observed in HTSC's are interpreted as arising from the 
sign change of the pair potential (PP) which leads to 
the formation || of midgap surface states. The spectral 
weight of these states has, for a d x i_ y i state, a maxi- 
mum for (110) oriented surfaces and vanishes for (100) 
surfaces. The absence of ZBCP in Ndi.ssCeo.isCuCU is 
considered evidence fl(]| for a dominant s-wave compo- 
nent. 

Previously published work on tunneling in unconven- 
tional superconductors examined the unpolarized case 
or |l|,|,| F/S junctions with an s-wave PP g§. 
Here we consider the theory of tunneling spectroscopy for 
an F/S junction with arbitrary spin polarization. We in- 
vestigate the interplay of ferromagnetism and unconven- 
tional superconductivity in forming ZBCP, ZBCD, and 
other features at finite bias. We include the effects of 
the exchange energy (related to the degree of polariza- 
tion) interfacial barrier height and Fermi wavevector mis- 
match (FWM) in the magnitude of the Fermi wavevec- 
tors in the F and S regions. Variation of these parameters 
leads to rich behavior and novel features in the conduc- 
tance, which require careful interpretation. Thus, FWM 



and spin polarization can combine to yield a ZBCP in 
an s-wave superconductor, while if one neglects FWM 
iflPU^JT^I the effect of spin polarization invariably leads 
to suppression of AR. 

We solve the Bogoliubov-de Gennes (BdG) equations 
@,|jl|Jl|l for a ballistic F/S junction. We extend the 
usual one-body Hamiltonian approach of Ref. Q to in- 
clude: 1) scattering at the F/S interface, x = 0, modeled 
by a potential V(r) — H5(x), where H is the strength 
of the potential barrier, and 2) allow for FWM, i.e. 
Ep = H 2 k F /2m in the F region at x < ( Ep is the 

2,2 



spin averaged value, Ep = (h kp*/2m 



h 2 k Fl /2m)/2,) 
and E' F = h 2 k' F /2m in the S region at x > 0. We in- 
clude the exchange energy Q h(r) — ho<d(—x), (Q{x) is 
a step function) and the pair potential A(k',r) = 

A(k', r)Q(x). From the invariance of the Hamiltonian 
with respect to translations parallel to x — 0, the par- 
allel component of the wave vector is conserved at the 
junction 11 ljj. Then, the parallel component of the so- 
lutions of the BdG equations is a plane wave, and the 
problem reduces to a ID one. 

For an electron injected from the F side, with spin 
S =t) I) excitation energy e, and wavevector kg at an 
angle 9 from the interface normal, there are (without 
spin-flip scattering) four scattering processes [ [ll]Jl4| with 
different amplitudes. For specular reflection at the in- 
terface, these are: 1) Andreev reflection J8|,[l3],[l5| with 
amplitude as as a hole with spin, wavevector, and an- 
gle with the interface normal, 5* (opposite to S), k-^, 
and respectively. 2) Ordinary reflection with ampli- 
tude bg as an electron with variables S, —kg, —6. 3) 
Transmission with amplitude c$ as an electronlike quasi- 
particle (ELQ) with k'cf and 9' s . 4) Transmission (am- 
plitude ds) as a holelike quasiparticle (HLQ) defined by 

FS 



k' 5 , and — Q' s . Here we have used H[Tl) kg^ik 



and k'^ « k' F , as explained below. The ELQ and HLQ 
have different, spin dependent, wavevectors, and there- 
fore they feel different pair potentials As + and Ag_ , with 
As± — \As±\exp(i(j)s±)- The spin dependence (even 
without FWM) of wavevectors and PP's is a novel feature 
of F/S junctions. 

Conservation of ku s yields the analogue of Snell's law, 



tp 5 sin t 



k' F sin i 



S' 



S =|, |, and for k FS > k' F 



1 



there is a spin dependent angle of total reflection. For the angularly averaged (AA) conductance jfl), < Gs >, 
a PP of the d x 2_ y 2 form, and allowing for different an- as < Gs >— J n d9 cos9G(9) / J n d9cos9, where fis is 



gles a e (— 7r/2,7r/2) between the crystallographic a-axis 
and the interface normal, we have A$± — Aq cos(29' s± ), 
where 9' s± are related to 9' s (and thus to the incident an- 
gle 9, from "Snell's law") by 0' s± =6' s Ta. The spin de- 
pendence of these angles will produce more complicated 
conductance features. 

In solving the BdG equations in the direction normal to 
the interface, we have for the magnitude of the relevant 
wavevectors: in the F region, fc^ = (2m/fi 2 ) 1 / 2 [Ep ± 
e + psho] 1 / 2 , where ps = ±1 for S =t (I), and in the 



S region, k 



2m/ti z ) 1 / 2 [E' F ± (e 2 



|2U/2jl/2_ 



limited by the angle of total reflection or by experimental 
setup. 

We concentrate on the d x 2_ y 2 state and take param- 
eter values appropriate for HTSC's and ongoing experi- 
ments jlTj on the effect of spin polarization on G. Typ- 
ically there is a small interface resistance || i.e. small 
Zq values, away from the tunneling limit Zq 3> 1 0. 
The experiments use H half-metallic ferromagnets with 
X -> 1 and FWM values of L < 1 @. We present 
our results for G and (G) as a function of dimensionless 
energy E = eV/Ag, where V is the bias voltage. We take 



In the regime of interest, Ep,E' F ^> max(e, |Ag±|) and E F /A = 12.5, with Ep = E' F (Lq = 1) and Ep = AE' F 



fc'c 



lUlql we have k Q 



(2m/h 2 ) 1 / 2 [E F + p s h \ 



1/2 



k' F . From this and "Snell's law", the components 



of kg and k^ normal and parallel to the interface can 



be found. We write k s = (ks, fcys), and k s 



t± — 



in the F and S regions. 

The conductance spectrum is calculated via the 
Blonder-Tinkham-Klapwijk (BTK) method @ ex- 
tended to include unconventional superconductivity ]Tl| ] 
and net spin polarization. It is sufficient QJl^l to calcu- 
late as and 6s: 



as 



AtsLsT+e-^ 3 * 



Uss+Ug S _ - Vss-V s - s+ T + r-eWs--^+) 



(1) 



b Vss+u-ss- - Uss-v-s s+ r + r-e^-^ 

s U ss +U- §s _-Vss-V- s - s+ r + r_e i ^—<l> s+ ) [> 

where we introduce T± = (e — (e 2 — \As±\ 2 ) 1 ^ 2 ) /\As±\, 
Ls = Lq cos0' s / cos9, Lq = k F /kp, describing FWM, 
t s = (l + PsX) 1 ? 2 , % = (1 - p s Xf/ 2 cos%/ cos6», X = 
ho/Ep , which defines the degree of polarization, f%5 ± = 
% + ws±, V S s± = t s - ws±, w s ± = L s ± 2iZ, Z = 
Zq/cosO, and Zq = mH/hkp is the interfacial barrier 
parameter. The normalized differential conductance p"l| ] 
(in units of e 2 /h) is then 

G ee G T + G i = W + ^l«s| 2 ~ \bs\ 2 ) (3) 

where the probabilities of an incident electron with spin 
5, P s , satisfy §, P T /P X = tt+X)/(l-X). At X = we 
recover the results of Rcf . |ll|] . The ratio of wavevectors 
in Eq. (||) reflects that the incident electron and the AR 
hole belong to different spin bands. One can use the con- 
servation of probability current fl4|| to generalize the sum 
rule for the reflection coefficients in the case of subgap 
conductance (E < |Ag±|), for the unpolarized case. We 
get g |a s | 2 + |6 S | 2 = 1 and then G s can be expressed in 
terms of the AR amplitude only. For fc^s > k' F , one sees 
from Eq. ([|) and "Snell's law", that for \9\ greater than 
the angle of total reflection (as — 0) Gs = 0. We define 



(L = 1/2), to consider the influence of FWM. 

In Fig. |l], we show results for G(E) at 8 = 0, and 
a = (i.e. an F/S interface along the (100) plane). This 
behavior is the same as for an s-wave superconductor 
(As± = Ao). In panel (a), at zero interfacial barrier, we 
display the effect of increasing X. The solid lines repre- 
sent results with FWM, L = 1/2. The behavior of the 
amplitude of G at zero bias (AZB) reflects the interplay 
between the effects of FWM and ferromagnetism. At 
X = there is a ZBCD, as in previous work [l^j . This 
is caused by the effective barrier introduced by FWM 
(even at Zq = 0) , at the interface which separates regions 
with different Fermi energies. With increased exchange 
energy the ZBCD evolves into a ZBCP, which narrows 
with decreasing Lq. The ZBCP in this case is not due to 
unconventional superconductivity, since there is no sign 
change in the PP's experienced by ELQ and HLQ. For 
reasonable values of the FWM, the maximum AZB is 2, 
independent of X. The AZB maximum is obtained from 
the condition k-\ki = kg = k' F . For Lq = 1/2, this occurs 
at X ~ 0.968. Thus, in the presence of FWM Andreev 
reflection can be enhanced by spin polarization and can 
even become maximal at a special value of X. These 
results differ from those obtained without FWM when 
AR is always suppressed and G decreases with X. The 
dashed lines (no FWM) illustrate this point. In panel 
(b) we show the influence of barrier strength at fixed 
X = 0.6. In the presence of FWM the subgap conduc- 
tance is more reduced, with sharper peaks at E = 1 than 
for the same Zq at Lq = 1 (dashed lines), since the effec- 
tive barrier strength Jfq] is enhanced. 

In Fig. U we use the same parameter values and no- 
tation as in Fig. |l| to display the effect of AA on G for 
s-wave (As± = Aq for all angles). In panel (a) the solid 
curves, (with FWM) show that the ZBCP at fixed X re- 
mains after AA, while its amplitude is typically reduced. 
We see that, unlike in the unpolarized case |18|, FWM 
can actually enhance G(0) at fixed polarization. In panel 
(b) we show that the effects of interfacial barrier on (G) 
are similar to those given in Fig. ^. 

We now turn to the effects of the sign change of the PP, 
i.e., of the unconventional nature of the S region. We take 
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a ^ and 6 ^ so that ELQ and HLQ may feel PP's of 
opposite sign. First, in Fig. [|, we consider the limit of 
no FWM, with 6 — it/ '12 and a = 7r/4. We examine the 
dependence of the AZB on X and Zo. In panel (a) we see 
that the ZBCD becomes more pronounced with higher X, 
because of suppressed AR. In the limit (Ep — ho) — » 0, 
the subgap conductance vanishes, due to the vanishing 
of the minority spin density of states (see the bottom- 
most curve). For X = 0.4 there appears a finite bias 
peak (FBCP), which moves to lower energy with increas- 
ing X. This follows from "Snell's law". At larger X, 
there will be increased difference between k$ and k' s , 
and 0' s will depart more from 8. ELQ and HLQ feel 
PP's with increasingly different spin dependent magni- 
tude: Gj- and Gj are governed by different energy scales. 
When a ^ 0, 7r/4, G(E) displays two or four distinct fea- 
tures determined by up to four different PP values. In 
panel (b) we show the decomposition of the X = 0.6 re- 
sult into its Gj and Gj components. The position of the 
FBCP, discussed above, is at the maximum of Gj. We 
also examine the effect of Zq at constant X. With in- 
creased Zq, the FBCP evolves towards smaller energies. 
Eventually, the barrier effects dominate those of X and 
the ZBCP resembles that found in the N/S junctions, 
attributed to midgap surface states ||[ll| . 

In Fig. [| we illustrate the effect of AA using the pa- 
rameters from Fig. ||. The solid curves represent AA over 
all angles below total reflection and the dashed curves are 
averages over a narrower region. Panel (a) displays the 
conductance for different polarizations and Zq — 0. The 
ZBCD reported in Q resembles the bottommost solid 
curve. The parameter values for this curve agree with 
their values in this experiment, as mentioned above. We 
see that the ZBCD in (G) occurs only at large X. The 
curves are then qualitatively different from those found in 
the tunneling limit for an s-wave superconductor, where 
the peak in (G) is sharp and at the gap energy. In panel 
(b) we show that formation of ZBCP with increasing Zq 
is a robust feature present in both types of AA. 

In Fig. H we consider the interplay of FWM and un- 
conventional superconductivity. We take Lq = 1/2. In 
panel (a) we show the results at Zq = for a range of 
values of X. The ZBCP evolves into a ZBCD. In panel 
(b) at X = 0.7, we see that with increasing Zq a FBCP 
forms and then evolves to a ZBCP. The FBCP here has 
a different origin than the breaking of the time reversal 
symmetry state in the N/S system [B 11 lj|. From Eqs. 

it is simple and instructive to obtain other G(E) re- 
sults. 

We have shown here that spin polarized tunneling spec- 
troscopy of F/S junctions displays qualitatively novel be- 
havior. The variety of features in G(E) and (G(E)) arises 
from the interplay among the form of the pair potential, 
the exchange energy of a ferromagnet modifying the AR, 
and the Fermi energies and interface properties of the F 
and S regions. The results are quite sensitive to FWM 



which should be carefully taken into account in inter- 
preting spin polarized experiments and can not, in the 
polarized case, be simply replaced by a change in param- 
eter Zq. 
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FIG. 1. G{E) (Eq. §). E = eV7A . Results are for 
8 = (normal incidence) and a = 0. The solid curves are 
for L = 1/2 (FWM present): in panel (a) at Z = (no 
barrier) they are (from top to bottom at E > 1) for exchange 
energies X = h /E F = 0,0.6,0.8,0.968,0.99,0.9999, while 
in panel (b) they are at X = 0.6 and (from top to bottom) 
Zo = 0, 0.25, 0.5, 1, 1.5. The dashed curves are for Lq = 1 (no 
FWM). In panel (a) X = 0,0.6,0.968, in panel (b) they are 
at Z = 0, 0.05. 



FIG. 2. (G(E)), the 6 averaged conductivity, for the same 
parameter values and curve identifications as in Fig. |l]. Re- 
sults in this Figure are averaged over all angles. 
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FIG. 3. G(E) for 6 = tt/12, a = tt/4, and 
I/O = 1. In (a), at Zo — 0, the curves are for 
X = 0,0.4,0.6,0.7,0.8,0.85,0.9,0.99, (top to bottom at 
E — 0.) In (b), the solid curves correspond, from top to 
bottom at E = 2, to Z = 0, 0.25, 0.5, 1, 1.5. The dashed and 
dash-dotted curves are, respectively, the Gf and Gj. conduc- 
tances at Zq = 0. 

FIG. 4. (G(E)) for the data in Fig. |[ The solid curves 
are averages over all 6 while the dashed curves are over a 
region of width 7r/24 centered at 7r/12. The curves shown 
are, respectively, in the same order as those in Fig. ^| except 
that in (a) results for X = 0.4, 0.7, 0.85 have been excluded 
for clarity. 



FIG. 5. G{E) for 6 = tt/12, a = tt/4, and FWM 
with Lo = 1/2. In (a), Zo = 0, and the curves are for 
X = 0,0.4,0.6,0.7,0.8,0.85,0.9 (top to bottom at E = 0). 
In (b), the influence of Zo is shown for X — 0.7. The 
curves correspond, from top to bottom at E = 2, to 
Z = 0,0.25,0.5,1,1.5 
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Fig 1, Zutic & Vails "Spin polarized.." 
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Fig 2, Zutic & Vails "Spin Polarized.." 
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Fig 3, Zutic & Vails "Spin Polarized.." 




Fig 4, Zutic & Vails "Spin Polarized.." 




